Strong-field approximation for Coulomb explosion of by short intense laser pulses 
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We present a simple quantum mechanical model to describe Coulomb explosion of H^" by short, 
intense, infrared laser pulses. The model is based on the length gauge version of the molecular strong- 
field approximation and is valid for pulses shorter than 50 fs where the process of dissociation prior 
to ionization is negligible. The results are compared with recent experimental results for the proton 
energy spectrum [I. Ben-Itzhak et al., Phys. Rev. Lett. 95, 073002 (2005), B. D. Esry et al., Phys. 
Rev. Lett. 97, 013003 (2006)]. The predictions of the model reproduce the profile of the spectrum 
although the peak energy is slightly lower than the observations. For comparison, we also present 
results obtained by two different tunneling models for this process. 

PACS numbers: 32.80.Rm,33.80.Rv,42.50.Hz. 
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I. INTRODUCTION 

Ultrashort, highly-intense, laser pulses at infrared 
wavelengths are currently being used to study molecu- 
lar dynamics under extreme conditions. These systems 
are commonly based upon a Ti:Sapphire laser, tunable 
around 790 nm, which can reach peak intensities in the 
range 10 13 -10 17 W/cm 2 , with pulse lengths shorter than 
50 fs. Consequently energy can be deposited on time- 
scales shorter than the fastest molecular vibration with 
field strengths comparable to the molecular bond. When 
molecules are exposed to such an environment, the re- 
sponse is highly-nonlinear and generally leads to mul- 
tiple dissociative ionization and high-order optical scat- 
tering. Indeed the multiple fragmentation process for a 
heavy polyatomic multi-electron system prevents a de- 
tailed analysis of the energy transfer, simply due to the 
proliferation of fragments produced. Small molecular 
systems, on the other hand, have simpler structure and 
fewer relaxation channels. Moreover, the fundamental 
molecules such as H2 or D2, have an intrinsic value owing 
to their fast vibration. Since a 790 nm pulse has a cycle 
period of 2.6 fs, and perhaps 10-50 fs duration, the vibra- 
tion provides an additional internal clock that records the 
response of the electronic excitation during the passage of 
the pulse. Consequently, there has been extensive study 
of the interaction of intense field dissociative ionization 
of H2 an with progressively shorter and more intense 
pulses and, in particular, using the ejected proton energy 
spectrum as an indicator of the electron response and a 
diagnostic of the pulse itself. For reviews of progress in 
this field one can consult, for example, Refs. [HELEl- 

The effect of the driving pulse on the rates of dissocia- 
tion and ionization as well as the energy of the fragments 
has been discussed in Refs. It has been noted that 

the kinetic energy release (KER) spectra indicated a de- 
crease in energy of the ionization products when the pulse 
duration was increased pi] . For a, comparatively long, 90 
fs pulse, the observation of a proton energy peak near 2 



eV would correspond to an internuclear separation R ~ 7 
ao at the instant of ionization. Given the equilibrium 
bond length of IFj" is R — 2 ao, this suggests dissociation 
occurring prior to ionization. The increased probabil- 
ity of ionization with increasing R can be understood by 
the mechanism of resonantly-enhanced multiphoton ion- 
ization via the antibonding 2pcr+ state 0, H| when the 
resonances are tuned 'from the red' as the bond expands. 
Conversely, at very high intensities the ionization process 
can be characterized as a tunneling transition. The as- 
sociated shape resonances of the molecular ion produce 
the effects known as charge-resonance enhanced ioniza- 
tion Q and dynamic tunneling ionization [Io| . These 
ionization mechanisms, diabatic (multiphoton) and adi- 
abatic (tunneling) ionization, can be studied indirectly 
via the proton emission spectrum. 

Perhaps the most interesting aspect of current experi- 
ments is that the laser intensities are high enough to ex- 
plore the adiabatic and diabatic regimes, and also pulse 
durations are short enough to probe the vibrational mo- 
tion. The study of bond-stretching during the ioniza- 
tion process has been made possible by enhanced con- 
trol over pulse duration [4j, |ll( and through pump-probe 
techniques [H EE 0. 

In Ref. [l3j it was shown that 
if H 2 was subjected to two ultra-short pulses, the rela- 
tive yield of low-energy protons raised at a probe delay 
of 24 fs. Similar evidence for dissociation prior to ion- 
ization was reported in Ref. 7] using single 100 fs pulses 
on a target. In this case, the KER-spectrum was 
explained by one- or two-photon absorption to the disso- 
ciation channel followed by ionization at a separation of 
12 a giving a structure reflecting the initial distribution 
of vibrational states. All of these experiments resulted in 
a maximum kinetic energy of the protons in the order of 
4 eV. 

To minimize the effect of dissociation and thereby be 
able to observe protons with higher energies, shorter 
pulses are needed. Experiments using 45 fs pulses showed 
KER-spectra in the range of 2-10 eV when the intensity 
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exceeded 1.7xl0 14 W/cm 2 @, E3, At lower intensi- 
ties the dissociation channel dominated completely and 
no Coulomb explosion was observed. The structure of 
these spectra was very recently qualitatively explained 
by associating the critical bond lengths with the Flo- 
quet (dressed) state [HI crossing with the ionic (1/R) 
curve. The underlying physical reason is that the den- 
sity of states is very high at threshold and there will be 
a resonant coupling to a quasi-continuum. The strength 
and width of the peaks are partly restricted by the num- 
ber of photons involved, and whether or not the pulse is 
long enough for the molecular expansion to occur, and 
so reach the crossing before the pulse ends. Using this 
model, and fitting a 12-parameter function, gives very 
good agreement with the observations. 

In this work we present a simple quantum mechani- 
cal model designed to predict the high-energy part of the 
KER-spectrum of dissociative ionization of . We use 
the strong-field approximation (SFA) within the length 
gauge O [13] • That is we calculate transition rates 
from a molecular state, in which the laser field is ne- 
glected, to a state in which the laser field is accounted 
for to all orders, but some approximations are made for 
the three-body continuum. The calculations are simpli- 
fied by assuming the Born-Oppenheimer separation of 
motions, and the Franck-Condon principle is invoked for 
the nuclear motion. These approximations allow us to 
determine the dissociative ionization rates to different 
channels in a simple way. The pulse is presented by a 
Fourier expansion so that the temporal intensity varia- 
tion is taken into account. This averaging is important 
due to the ponderomotive shift of the electron which may 
cause channel closings when the intensity is raised. 

Following this procedure, and integrating over the elec- 
tron spectrum, the proton energy distribution can be ob- 
tained. Comparing with recent experiments pi [l5|. we 
find the model predictions to be in very good qualitative 
agreement for the shape of the distribution. However, our 
calculations predict the peak of the proton spectrum at 
slightly lower energies than that observed in experiment. 
Finally, for completion, we compare the predictions of 
our model with tunneling theory calculations. Atomic 
units (|e| = h = m e = ao = 1) are used throughout 
unless indicated otherwise. 



II. MODEL 

We consider the transition from an initial field-free 
state to a state in which the electron is only affected 
by the laser light, and the two protons are only subject 
to their mutual Coulomb repulsion. There are three es- 
sential elements of the model. Firstly, the initial state, 
in which the vibrational population distribution plays 
a key role. Secondly, the overlap of these vibrational 
states with the final Coulomb states of the proton pair 
is of great importance in modulating the proton spec- 
trum. Thirdly, and most importantly, the coupling to 



the continuum, the nonpertubative photoionization rate, 
as a function of bond length and laser intensity will de- 
termine the range of proton energies. All three factors 
are intrinsically linked. 

Consider a monochromatic plane-wave component of 
the light field, with linear polarization, z, and angular 
frequency uj so that the vector potential in the dipole- 
approximation is 

A(t) = A zf(t)cosujt , (I) 

where Ao is the field amplitude, and the pulse shape is 
described by the factor, < f(t) < 1, which we take to 
be a Gaussian profile. 

A. Molecular states 

The process of formation of , requires ionization of 
the neutral species. The molecular ion H^" has only a sin- 
gle bound electronic state (ls<7 fl ) that is exactly known. 
The nuclear relaxation that follows this primary ioniza- 
tion is not so well defined and remains a source of uncer- 
tainty and investigation |2j. Nonetheless, to a very good 
approximation, the rotational degrees of freedom can be 
considered as frozen since its characteristic timescale of 
~ 170 fs is much longer than the pulse duration. How- 
ever, this results in an ensemble of vibrational modes, 
as observed in experiment [l8j| . As usual, the z-axis of 
the laboratory reference frame is defined by the polariza- 
tion vector, while the z-axis of the molecular (body- fixed) 
frame is defined by the internuclear axis. The notation 
for the coordinates is that R denotes the internuclear vec- 
tor and the electron coordinate with respect to the inter- 
nuclear midpoint, is denoted by r. In the following pre- 
sentation, for consistency, the laboratory reference frame 
is employed. If we let v denote the vibrational quantum 
number, then the eigenstates of the initial ensemble can 
be expressed as: 

r, = Mr, R)xUR)e- iEi "\ (2) 

where (f>i(r,R) is the electronic wave function, with en- 
ergy £i(R), and Xiu{R) is the vibrational cigenfunction 
with eigenvalue E^w- The total energy is then, Ei V — 
£i(i?o) + Enw- Recall that the electronic function has a 
(7+ symmetry, and the transformation of 4>i(r,R) from 
laboratory to molecule frame is effected by the Wigner 
D rotation matrix [l9j . 

In the final state, given the large separation between 
the nuclei and the electron at the time of ionization, we 
suppose that the influence of the protons on the electron 
is negligible compared to that of the external field [TtJ • 
This is consistent with the asymptotic (t — > +oo) limit of 
the system as a (decoupled) product state of an outgoing 
Volkov wave 4>f (electron in the electromagnetic field) 
and a Coulomb wave Xf f° r the proton motion: 

* / (iJ,r > t) = / (r,*)x/(ie,«). (3) 
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For an electron in a laser field described by (JTJ) , the 
length-gauge Hamiltonian is given by 



m 



r-F, 



(4) 



where p is the canonical momentum and, F = —dtA, 
represents the electric field. The Volkov states form a 
complete set of solutions to the equation Hj lec (f>f(r, t) — 
idt<pf(r,t) and can be written: 



<t> f {r,t) = exp 



i(q+ A (m . r -if <i±4«>V 

J -oo * 



B. Transition rates 

To derive the expression for the transition amplitude, 
we follow the procedure of Ref. 20], and generalize to 
the case of an incoherent mixture of initial vibrational 
states each populated with probability P v , defined by 
the Franck-Condon factors. In this way we obtain the 
rate w for a transition into the final state ^ v of Eq. ([3]). 

oo 



,(5) 



where q is the kinematic momentum corresponding to 
a drift energy q 2 /2, and we denote the ponderomotive 
energy as U p = Aq/4. 

Since vibrational energies are much larger than rota- 
tional energies, and Coriolis coupling can be neglected 
at these energies, we make the usual assumption that 
the proton ejection occurs along the internuclear axis 
direction without rotation: the axial recoil approxima- 
tion. Thus, the nuclear motion is governed by the one- 
dimensional Coulomb repulsion: 
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where fi = \m p is the reduced mass. The corresponding 
eigenfunction, with energy -EV/, and wavenumber, kf = 
has the form, 



(7) 



where 



■ exp 



7T fl 
2 Vf 



+ ik f R 



Til + i 
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k f R 
(8) 



XxFx(l + ij^;2;-2ik f R), 

and \F\ is the confluent hypergeometric series. Since the 
Volkov wave represents all orders of the field amplitude, 
the final state is a coherent sum of the full spectrum 
of photoelectron harmonics, including the angular distri- 
bution. Since our primary interest here is the compari- 
son with experiment for the proton energy spectrum, we 
do not present results for the photoelectron differential 
yields. Instead we must integrate over these degrees of 
freedom. Furthermore, the initial vibrational state is a 
mixed state and the orientation of the molecule is ran- 
dom. In all, this amounts to integrating over four contin- 
uous variables and summing over two discrete variables, 
in addition to the matrix clement calculation. However, 
since numerical quadrature is inherently a set of inde- 
pendent calculations, these calculations can readily be 
performed on a parallel computer. 
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where 



Aun=^ f (^ f \r-F\^ iv )dt, 
1 Jo 



(10) 



and the minimum number of photons absorbed, no, is de- 
termined by energy conservation. For a given number of 
absorbed photons, n, the electron momentum is defined 
as: 



q = J 2{E iv +ntu -Up- E Nf ). 



(11) 



To obtain the dissociative ionization rate, we multiply 
by the density of states per unit energy, per unit solid 
angle. Using the normalization convention defined above, 
the appropriate factor is (2ir)~ 3 qd 2 q. So that we have: 



dw 
dE N f 



i'7l\ J 



(12) 



where dq defines the direction of the outgoing electron. 

The calculation of A un {q) is significantly simplified in 
the Franck-Condon approximation, where it is assumed 
that the electronic transition appears almost instanta- 
neously compared to changes in the nuclear position. 
That is, we can make the integration over the electron 
coordinate independent of the nuclear coordinate by re- 
placing R by some fixed value Rq, 

i r T 

A m =S fi - D e ,(Ro,t)e* E »t- E "'»dt, (13) 
1 Jo 

where the Franck-Condon factor S/i and the electronic 
matrix element D e i (Rq , t) is given by 

S fi = f X *f(R)xMdR, (14) 
Jo 

D el (Ro,t)= J cj ) }(r,t)r-Fcj )l (r,R )d 3 r. (15) 

We have tested the validity of the Franck-Condon ap- 
proximation and found it to be accurate, especially for 
small v. 

The evaluation of the matrix element is conveniently 
carried out in spherical coordinates [l7| • The first step is 
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to rewrite the expression as 



D el (Ro,t) = E lv -E_ 



<J Nf 



(q + A) 



x exp 

where we have used 



dt 



r ■ F 



(16) 
iu(q + A,Ro), 

(17) 



the axial-recoil model, to summing over the ejected pro- 
ton directions described by the rotation matrices. This 
finally reduces to an energy-differential electronic rate, 
T v (E N f), 



and 4>i denotes the Fourier transform of the electronic 
wave function, 



i{k,Ro) 



exp [—ik ■ r] 4>i(r, Ro)d 3 r. 



(18) 



Here k = q + A. In the length gauge formulation of 
the SFA, the transition amplitude only depends on the 
asymptotic form of the coordinate space initial electronic 
state [ljl [13, H3, Hl| • In the laboratory fixed frame this 
electronic wave function for nuclear orientation Rq reads 



i(r,Ro) 



- 1 ) exp (-«r) ]T £ C l0 Dl2 Q (Ro)Y lm (r), 



I m 



(19) 



with Ciq asymptotic expansion coefficients [22| . and 




(20) 



where is the eigenvalue of the electronic Hamiltonian 
including the nuclear repulsion. In Eq. (fl~9|) D^ (Rq) is 
the Wigner rotation function that effectuates the trans- 
formation from the molecular to the laboratory fixed 
frame. 

Finally the electronic matrix element can be found as 



D el (Ro,t) = E,, 



E Nf - 



(q + A) 



x exp 



2 



l m 



oo 



: / ji(kr)r^ K ^ exp (— nr) r 2 dr. (21) 
/o 

The radial integral has a closed analytic form in terms 
of Gauss's hypergeometric function [23[. The time- 
integration is performed numerically, along with the sum 
over the number of photons. In the experimental spec- 
tra, the protons resulting from dissociative ionization are 
collected over a range of ejection angles, with respect to 
the polarization direction. Averaging over the different 
orientations of the molecular axis is equivalent, within 
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D el (R ,t)e l( - EN f' E ^ t dt 



(22) 



To obtain the total rate we use Eqs. p^ - ([T3)) . i.e., we 
multiply by the Franck Condon factor 



dw y 



\Sfiv(E N f)\ 2 r„(E N f) 



dE Nf 

and sum over the different initial states 

dw dw v 
dE Nf ~ 2-*> v dE Nf ' 



(23) 



(24) 



The values of P v are determined by the formation mech- 
anism of the ion. In the context of recent experiments 
[3, [lH it is reasonable to assume this to be a Franck- 
Condon distribution. 

Finally, to compare with experimental data, we aver- 
age over the pulse profile, f(t) [Eq. (JT])], which is taken 
to have a Gaussian profile with FWHM 45 fs. Under 
the assumption that the variation in the pule envelope 
is slow compared with the optical period, the definition 
of ionization rate for fixed Aq is still valid. The ion- 
ization process can be significant for intense pulses, and 
thus we should allow for depletion of the molecular state. 
To model this process, the total ionization probability is 
found by integrating this rate over time 



Pj(E Nf )=J2 



dw v (t) 
~dE 



Nf 



N v {t)dt. 



(25) 



Here N v {t) denotes the population in a given vibrational 
state v. This population can be found from the rate 
equation: 



dNjf) 
dt 



= -T» ot N v {t), 



(26) 



Pv and Y\ ot 



with the boundary condition N„(—oo) 
denoting the total rate of ionization from the vibrational 
state v. By integrating the differential rate over all final 
states this rate is found as 



dw v {t) 
V tot~ I — dE Nf- 



dE 



Nf 



III. RESULTS 
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In panels (a) and (c) of Fig. Q] we show the predicted 
KER-spectra for Coulomb explosion of H^ for two pulses 
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FIG. 1: The predicted KER-spectrum for pulses character- 
ized by a wavelength of 790 nm, a duration of 45 fs (FWHM) 
and a peak intensity of 4.1 x 10 14 W/cm 2 (a) and 8.7 x 10 14 
W/cm 2 (c). We also plot the experimental results obtained 
with these pulses in ji| (panel (a)) and [l5| (panel (c)), re- 
spectively. Panels (b) and (d) show the contributions to the 
overall spectrum from the different initial vibrational states. 
The darker the shading the larger the contribution to a given 
final proton energy from a particular initial vibrational level. 



used in recent experiments 0, [l5j . The spectrum corre- 
sponding to a pulse duration of 45 fs and a maximum 
intensity of 4.1 x 10 14 W/cm 2 peaks near 3 eV, while the 
spectrum is moved to slightly higher energies for pulses 
with maximum intensity of 8.7 x 10 14 W/cm 2 . Along with 
the theoretical predictions, experimental results are also 
given [j, EH • The qualitative shape of the distributions 
is in good agreement in Fig. [T] (c). However, the exper- 
imental results at the lower intensity, Fig. [T] (a), show 
a much broader distribution, with a significant amount 
of fast protons (energies above 6 eV) not predicted by 
theory. 

To explain the structure of the theoretical spectra we 
study the contributions from the different initial vibra- 
tional states. Since each component contributes inco- 
herently, they can be studied in isolation, and we plot 
in Figs. [Ub) and (d) these contributions to the overall 
spectrum. Here we see that even though the molecu- 
lar ions are predominantly in low vibrational states, the 
bulk of the Coulomb explosion yield is from highly vibra- 
tionally excited states. For pulses of peak intensities near 
4.1 x 10 14 W/cm 2 ionization from the vibrational states 
v = 6— 11 is favored, while the v — 4, 5 and 6 states dom- 
inate for pulses of peak intensities near 8.7 x 10 14 W/cm 2 . 
It is this shift in origin of the protons that causes the shift 
in energy of the proton spectrum. Highly vibrationally 
excited molecular ions contribute to the low energy part 
of the spectrum, while low-excited molecular ions result 
in protons of a higher energy. 

To see how the structure of the contributions comes 
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FIG. 2: Calculations showing the result for ions in a vibra- 
tional excited initial state with v = 2. The initial and final 
wave functions squared are shown in panels (a) and (b). The 
energies written in panel (a) indicate the proton energy that 
the particular states correspond to (i.e., one half of the energy 
of the state). Panel (c) shows the Franck-Condon overlap be- 
tween the initial vibrational state and the final Coulomb wave 
for the internuclear coordinate. Panel (d) shows the elec- 
tronic ionization rate for two different intensities and panel 
(e) shows the ionization probability using pulses of 45 fs du- 
ration (FWHM), 790 nm wavelength and a peak intensity of 
8.7xl0 14 W/cm 2 , again for v=2. 



about, we examine the initial and final nuclear wave 
functions. In panel (a) of Fig. 2 the probability densi- 
ties of two energy-normalized continuum states are pre- 
sented, along with the density of the v = 2 state in 
panel (b) . The well known structure of v + 1 peaks with 
large probabilities near the classical turning points is seen 
here. Since we have a continuum of final states which 
are all peaking at different internuclear separations, this 
initial structure is reflected in the Franck-Condon over- 
lap Sfi V {EMf) shown in Fig. 2(c). Remember that large 
internuclear separations R correspond to low energies. 
In Fig. 2(d) we show the electronic transition rates for 
two different intensities. The rates decay, roughly ex- 
ponentially, with increasing ejected proton energy with 
an attenuation more pronounced for the higher inten- 
sity. When multiplying the Franck-Condon factors and 
the electronic transition rates and integrating over time 
we obtain the final result for the v = 2 case in Fig. 2(e). 
The influence of the Franck-Condon factor is still appar- 
ent but the attenuation produced by the decrease in elec- 
tronic transition rate suppresses all but the low energies 
in the proton spectrum. This is reflected in Figs. 1(b) and 
1(d), which isolate the contribution of the different vibra- 
tional states weighted by their populations. The domi- 
nant peak in the contribution is moved towards lower 
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energies as the molecular ion gets vibrational excited, re- 
flecting the structure of the nuclear wave function show- 
ing a classical turning point at a larger nuclear separa- 
tion. 

Comparing Fig. 1(b) and (d) shows that at the higher 
intensity each vibrational component is broader and the 
combined spectrum extends to higher energies, which 
is also expected due to the less rapid decrease of the 
electronic ionization rate at this intensity. In addition 
the largest contributions now come from less excited vi- 
brational states. The combination of increase in the 
high-order multi-photon ionization rate and the relatively 
larger population in the lower v states compared to the 
higher v states leads to this enhancement of faster pro- 
tons. 

Overall we see a movement of the predicted spectrum 
towards higher energies as the intensity is raised, both 
due to the change in the contributions from the differ- 
ent vibrational states and due to the favoring of low ex- 
cited vibrational states. This tendency is reflected in the 
measurements, but the resemblance is far from perfect. 
The lack of agreement may be attributed to several fac- 
tors. Clearly, a limitation of our model is the calculation 
of ionization using the strong-field model which is rele- 
vant to adiabatic quasi-tunneling. If resonant-enhanced 
ionization at smaller bond lengths were significant, this 
would produce faster fragment protons. Also rescatter- 
ing effects are ignored, which might as well give protons 
of high energy. 



A. Results using tunneling 

To complete the discussion of Coulomb explosion of 
H^" the results using a simple tunneling theory are given. 
Here the frequency dependence of the field is neglected, 
and the laser is treated as a slowly-varying (quasi-static) 



electric field that produces tunneling ionization. Strictly 
speaking, this static model is only valid when the tun- 
neling time is much shorter than the optical cycle as ex- 
pressed by 7k *C 1, where jk — *J ^ s the Kddysh 

parameter. For I = 4.1 x 10 14 W cm -2 , jk ~ 0.8, while 
for I = 8.7 x 10 14 W cm" 2 , lK ~ 0.5. 

In the tunneling limit the ionization rate at a specific 
nuclear separation is given as [22| 



3F o y B 2 (m) 1 



,.3 



x exp 



ttk 3 J 2M|m|! k -Bt V 
2k 3 



3Fr 



»Tl|-l 



(28) 



where Fq is the electric field amplitude, Z = 2 for H 
K = y/2I p (R) and 



2 • 



B(m) 



-J2ci D l m0 (R)(-iy 



' (21 + l)(l + \m\)\ 



2(1 - \m\)\ 



(29) 



According to [22j the effect of nuclear motion is included 
by weighting the electronic ionization rate at different 
internuclear separations with the probability of being at 
this separation (the reflection principle); 



dw 
dR 



r tun (R)\ X i(R)\ 



(30) 



This result can be translated into a function of the proton 
kinetic energy by assuming that all the Coulomb energy 
between the two protons at the time of ionization is trans- 
lated into proton ejection energy. The results are shown 
in Fig. 3. 

Compared to the strong-field approximation model the 
low energy part of the spectrum shows better agree- 
ment with experiments, but the high-energy part is again 
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poorly represented. The reason for the better agree- 
ment in the low-energy part of the spectrum is clear 
from Figs. 3(b) and (d): a shift towards lower vibrational 
states and an exponential favoring of the low energy part 
of the spectrum as is typical for a tunneling theory. The 
lack of agreement at high proton energy could be an indi- 
cation that multiphoton resonance ionization is involved 
at smaller bond lengths, and this is clearly absent within 
the static field tunneling model. 

Another way of using tunneling theory to describe 
Coulomb explosion of is to include the effect of nu- 
clear motion as described in Ref . [25| . The overall rate is 
here found by 



dw 
dE Nf 



Xf(R)Tl n (R) Xi (R)dR 



(31) 



The results are given in Fig. 4, and show an even larger 
shift towards lower vibrational excited states resulting 
in a slightly better agreement particularly for the high 
intensity case. 



IV. CONCLUSION 

In this work we have presented a relatively simple 
quantum mechanical model to describe Coulomb explo- 
sion of molecular hydrogen ions. The model builds on the 
length gauge molecular strong-field approximation and 
gives reasonable predictions when compared with exper- 



iment. Theory predicts a dominant proton emission at 
slightly lower energies than the experimental measure- 
ments. However, the effect of multiphoton resonance 
ionization contributions, known to be significant for this 
range of intensities, could be important. 

The strong-field approximation allows very simple and 
fast calculations and might, for that reason, be a useful 
tool in the further understanding of molecular dynamics, 
including related intense field processes as, e.g., high- 
harmonic generation. For all purposes it is important 
only to use the model in the regime, where it is valid, 
namely describing pulses of high intensities and short du- 
rations. If the pulse duration exceeds 50 fs the molecular 
ion will have time to dissociate and the resulting KER- 
spectrum will move to considerably lower energies. 

For comparison we have discussed the predictions from 
two different tunneling models. The spectrum is here 
somewhat narrower and also centered near to low ener- 
gies. 
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